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I. INTRODUCTION 



Alan Coley* and Martin Goliath* 

Closed, spatially homogeneous cosmological models with a perfect fluid and a scalar field with 
exponential potential are investigated, using dynamical systems methods. First, we consider the 
closed Friedmann-Robertson- Walker models, discussing the global dynamics in detail. Next, we 
investigate Kantowski-Sachs models, for which the future and past attractors are determined. The 
global asymptotic behaviour of both the Friedmann-Robertson- Walker and the Kantowski-Sachs 
models is that they either expand from an initial singularity, reach a maximum expansion and 
thereafter recollapse to a final singularity (for all values of the potential parameter k), or else they 
expand forever towards a flat power-law inflationary solution (when k 2 < 2). As an illustration of 
the intermediate dynamical behaviour of the Kantowski-Sachs models, we examine the cases of no 
barotropic fluid, and of a massless scalar field in detail. We also briefly discuss Bianchi type IX 
- models. 

o 
o 

(N 

Oh, Cosmological models with a scalar field and an exponential potential are of fundamental importance in the study of 
the universe. These models are motivated by the fact that they arise naturally in alternative theories of gravity and 
occur as the low-energy limit in supergravity theories . By using qualitative techniques, the well-known power-law 
inflationary solution has been shown to be an attractor for all initially expanding Bianchi models (except a subclass 
of the Bianchi type IX models which will recollapse) in the class of spatially homogeneous Bianchi models . More 
-H , recently, cosmological models which contain both barotropic matter and a scalar field with an exponential potential 
have been studied Q , partially motivated by the fact that there exist spatially flat isotropic scaling solutions in which 
the energy density due to the scalar field is proportional to the energy density of the perfect fluid |J. In JtJ the 
stability of these cosmological scaling solutions within the class of spatially homogeneous cosmological models with a 
perfect fluid subject to the equation of state p — (7 — l)/z (where 7 is a constant satisfying < 7 < 2) was studied and 
it was found that when 7 > 2/3, and particularly for realistic matter with 7 > 1, the scaling solutions are unstable; 
essentially they are unstable to curvature perturbations, although they are stable to shear perturbations. In addition, 
in j|] homogeneous and isotropic spacetimes with non-zero spatial curvature were studied. 

It is clearly of interest to study more general cosmological models. One class of models of particular interest are 
Q-( those with positive spatial curvature. These models have attracted less attention since they are more complicated 
mathematically. Positive-curvature Friedmann-Robertson- Walker (FRW) models P~p3|, Kantowski-Sachs models 
(^jqI jl2| and Bianchi type IX models |3 101^,|lJ] have been studied using qualitative methods, although rigorous analyses 
using a new set of compact variables have not been carried out. The Bianchi type IX models are known to have 
.J^ j very complicated dynamics, exhibiting the characteristics of chaos [ ^0|JT5| ], and are hence beyond the scope of the 
present study. Recently positive-curvature FRW models and Kantowski-Sachs models with a perfect fluid and a 
cosmological constant have been investigated using qualitative methods and utilizing compactified variables. 

The outline of this paper is as follows. We shall first comprehensively study the qualitative properties of the class 
of positive-curvature FRW models with a barotropic fluid and a non-interacting scalar field with an exponential 
potential, extending and generalising work by Turner fl7j| , who used different basic variables. We shall then analyse 
the qualitative properties of the Kantowski-Sachs models. Positive-curvature FRW models and Kantowski-Sachs 
models belong to the class of spherically symmetric models, and hence the present work is a natural extension of 
recent work |lj||l9| ]. Indeed, it turns out that understanding the dynamics of the Kantowski-Sachs models is crucial 
for understanding the global dynamics of general spherically symmetric similarity models p0|. 
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A. Matter model 



The matter content of the models is taken to be a perfect fluid and a scalar field with exponential potential. The 
corresponding energy-momentum tensor is 

Tab = (Tp{ ab + Tgfab), (1) 

Tplab = nu a u b + p (u a u b + g a b) , (2) 

T s iab = 4>,a4>,b - + V ^^j 9ab, (3) 

V(cf>) = V e- K t (4) 

where k is a non-negative constant, and the pressure is given by p = (7 — with the equation-of-state parameter 
in the range 1 < 7 < 2. The fluid energy density fj, and the scalar field <j) are functions of a timelike coordinate t. 
A dot denotes differentiation with respect to t and throughout units are used in which c = 8irG = 1. The matter 
components are assumed to be non-coupled, and thus they are separately conserved: 

V a T p f = 0- V a Tjf. (5) 

For convenience, we define 

X = -Lj>. (6) 

II. CLOSED FRIEDMANN MODELS 

We start our investigation of closed cosmological models with a perfect fluid and a scalar field by looking at the 
closed FRW models. The line element for these models can be written 

ds 2 = -dt 2 + S{t) 2 dr 2 + S(t) 2 sin 2 rdn 2 . (7) 

The expansion of the fluid congruence is given by 6 — 35/5, and the evolution equation for the curvature K = 9/(9 2 S 2 ) 
is 

K = - — (38 + 8 2 ). (8) 



The conservation equations yield 



fi = -70/z, (9) 

x = -ex + (10) 



From the field equations we obtain 



n= l -[{\ + K)e 2 -ix 2 -w], (11) 
§ = -\{ e2 + \° 2 + \ t 3x2 " w + 3(7 " } • (12) 



Assuming /1 > 0, the Friedmann equation, Eq. (|il|), shows that D = + K)0 2 is a dominant quantity. Thus, 
compact variables can be defined according to 



9 Vsx Tlr VW 

Qo = ^, U =nr , W = —. (13) 



Note also that the curvature is given by 



K= 1 -^. (14) 
^0 
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The Friedmann equation becomes 



n D = ^ = i-u 2 -w 2 . 

Defining a new independent variable, ' = d/dr — d/dt, the evolution equation for D 



(15) 



D' 



-3Qo ( u z + 



D 



decouples. Thus, a reduced set of evolution equations is obtained: 



(16) 



Qo = (i-Oo) [i - 3 (y 2 + ln D ) 



U' = 3Q U 



-1 



U 2 + 1si D 



W = 3Q (u 2 + 



w 



(17) 



^kUW. 



There is also an auxiliary evolution equation 

fi'r> - -3Qo [(1 ~ Vnh - 2U 2 ] Q D , 



(18) 



and it is straight-forward to consider the set of variables (Qo, U, £Id), rather than (Qo, U, W) [jl7|. Note that by setting 
K = 0, U = 0, and identifying A = Vq, J7a = W 2 , the evolution equations corresponding to closed FRW models with 
a cosmological constant are obtained (Hi) . 

It is also useful to consider the deceleration parameter, given by 



q P i = -\ l 



UpfVqflpf 



1 

Ql 



1-3 U 



7 



(19) 



for Qo ^ 0. From this expression, we can see that there is an inflationary region (q p f < 0) in the state space whenever 
Q.D < — 3U 2 ). However, as will be seen below, it is only when k 2 < 2 that there exist attractors that are 
inflationary. Note also that for Qo 7^ 0, 



Qo = -(l-Qo)Qo9 P f, 



(20) 



so that Qo < whenever q p f > in which case Qo is itself monotonic. When q p { < 0, that is in the inflationary 
region, Qo need not be monotonic - for example, see the orbits close to +<& in Figs. |] and y. 
The dynamical system Eqs. (O) is symmetric under the transformation 



(t,Qo,U,W)^(-t,-Q ,-U,W) 



(21) 



Thus, it is sufficient to discuss the behaviour in one part of the state space, the dynamics in the other part being 
obtained by Eq. (^l|). 
Furthermore, note that 



M 
M' 



-(3 7 -2)Q M 



(22) 



is a monotonic function in the regions Qo < and Qo > for tip ^= 0. As there are no equilibrium points with 
Qo = when 7 > 2/3, M acts as a monotonic function in the interior of the state space. Consequently there can be 
no periodic or recurrent orbits in the interior state space and global results can be deduced. In addition, from the 
expression for the monotonic function M we can see immediately that either Qq — > 1 or — > asymptotically. 
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Interpretation 


Qo 


u 


w 


Note 


±F 
±K± 


Flat Friedmann 
Kinetic dom. 


€ 


n 

u 

±1 


n 







Scalar-field dom. 


e 


— f 




k 2 < 6 


±x 


Curvature scaling 


K c 

75 e 


7! 




< 2 


±FS 


Flat matter scaling 


e 


,/Il e 

V 2 k 


^1^7(2-7) 


> 37 


s± 


Static 





I / 3 7 -2 
3(2-7) 





7 <2/3 



TABLE I. Equilibrium points of the closed FRW models. 









<7 P f 


±F 


I 





^(37-2) 


±K± 





I 




±$ 





I 


±( K 2 -2) 


±x 





I 





±FS 

s± 


^(« 2 -3 7 ) 

4 

3(2-7) 


37 
7* 
37-2 
3(2-7) 


1(37^2) 
Not defined 



TABLE II. The physical quantities ilo, fi^, and q p t for the different equilibrium points of the closed FRW models. 



A. Equilibrium points of the closed FRW dynamical system 

A number of equilibrium points can be found for the dynamical system, Eqs. (0). In what follows, e = ±1 denotes 
the sign of Qo, whereas VL^ — U 2 + W 2 is a density parameter associated with the scalar field. In Table |, the various 
equilibrium points are summarized. The subscripts on the labels have the following significance: The left subscript 
gives the sign of Qo and indicates whether the corresponding model is expanding (+) or contracting (— ). The right 
subscript gives the sign of U; i.e., the sign of (f>. 

The equilibrium points labeled ±F correspond to the flat Friedmann solution. For these points, the scalar field 
vanishes ({7 = = W). There is an orbit from + F to _F and this orbit represents the closed FRW solution with no 
scalar field, starting from a Big Bang at + F and recollapsing to a "Big Crunch" at _F. 

The K points represent exact solutions with a massless scalar field (W — 0). As the fluid is negligible (tip = 0), 
these solutions are dominated by the kinetic term U. They correspond to Jacobs analogues of Kasner solutions in 
which U takes on the role of a shearing mode pj. 

There are equilibrium points with non- vanishing potential, where the scalar field dominates (flo = 0). These points 
±$ are only physical when n 2 < 6. For k 2 — 6, +$ coincides with +K + , and _<E> with _K_. For k 2 > 6, they are 
outside the physical part of the state space. The equilibrium point with e = +1 (i.e., + $) is a sink, and for k 2 < 2 it 
corresponds to the power-law inflationary attractor solution (cf. the expression for q p { in Table 

There are also points ±X for which the matter is unimportant, but the curvature is non-vanis: 
tracks the scalar field. The corresponding solutions are called curvature scaling solutions ||. These solutions only 
exist when n 2 < 2. For k 2 = 2, + X coincides with + $, and _X with Above this value of k, these equilibrium 
points are outside the physical part of the state space. 

When k 2 > 37, there is a flat matter scaling solution, for which both the fluid and the scalar field are dynamically 
important. The corresponding equilibrium points are denoted ±FS, and for k 2 = 37 they coincide with ±$. 

Finally, there are equilibrium points S± corresponding to static solutions, analogous to the Einstein static universe. 
These are only physical when 7 < 2/3. For 7 = 2/3, a set of equilibrium points appears along the line U = = W, 
signaling a change of stability when the points S± leave the physical state space. In what follows, we will only consider 
equations of state for which 1 < 7 < 2, hence we will not consider these points further. 

Table |0| presents some physical quantities for the various equilibrium points, and Table III lists their eigenvalues. 
The sources and sinks of the dynamical system when 7 > 2/3 are listed in Table |v|, and the global behaviour for 
different values of k can be summarized as follows: The state space when < k 2 < 2/3 is depicted in Fig. [l], where 



ting (Qo 1), and 
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Eigenvalues 



3e — J ffisgn (U) 



±F 
±K± 
±* 
±X 

±FS 

S± 



(3 7 - 2)e 

-^(37-2)6 
(3 7 ^ 2)e 



-|(2 - 7 )e 
3(2 - 7 )e 
4(6-« 2 )e 



-(2-« 2 )e 



"75 V~ 



^sgn([7) 



-I [(2 - 7) ± i ^(2 - 7 )[24 7 2 - (9 7 - 2)«*]] e 
±\/2V-(37-2) 



TABLE III. Eigenvalues for the different equilibrium points of the closed FRW models. 



Past attractors 



Expanding from a singularity (Qo > 0) 


+K_ 


Always 






+K+ 


k 2 < 6 




Contracting from a dispersed state (Qo < 0) 


_<E> 


k 2 <2 (and k 2 


<3 7 ) 


Future attractors 


Contracting to a singularity (Qo < 0) 


-K+ 


Always 






_K_ 


n 2 < 6 




Expanding to a dispersed state (Qo > 0) 


+ <& 


k 2 <2 (and k 2 


<3 7 ) 



TABLE IV. Summary of sources and sinks for the closed FRW models. 
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FIG. 2. The state space for closed FRW models with a scalar field when 2/3 < k 2 < 2. Dashed curves and white arrows and 
circles are screened. Dotted orbits are in the interior of the state space. 




the features in the rear part of the state space have been suppressed. Dashed and full curves represent orbits in the 
boundary submanifolds, while dotted curves represent orbits in the interior. Both +K_ and +K + are past attractors, 
while _K + and +$ act as future attractors. There are also orbits from +K_, whose future attractor is _K_ at the 
rear of the figure. Note that orbits future asymptotic to + <I> correspond to solutions that exhibit power-law inflation 
(— 1 < q p { < 0, see Table ||). Observe that the outgoing eigenvector directions from the saddle point + F span a 
separatrix surface in the interior of the state space. Similarly, + X is a saddle for which the ingoing eigenvector 
directions span another separatrix surface. These separatices confine orbits in the interior state space to specific 
regions. For example, there is one region where all orbits are past asymptotic to +K + and future asymptotic to -K + . 

When k 2 > 2/3, the separatrix surface associated with + X changes structure, see Fig. |^. This is easiest seen by 
considering the separatrix orbits in the flu — submanifold, see Fig. [| When k 2 < 2/3, there is one separatrix orbit 
from to + X, and one from _X to +$. For k 2 = 2/3, these two orbits coalesce into a single orbit from _X to + X. 
This orbit corresponds to a special "Bouncing Universe" solution, existing only for this particular value of k. When 
k 2 > 2/3, the separatrix orbit to +X starts at +K_, and the orbit from _X goes to -K + . Thus, there is a bifurcation 
for k 2 = 2/3; however, we note that there is no stability change of equilibrium points involved. A similar behaviour 



of separatrix surfaces has been found for Bianchi type IX models |13 



When k increases, the equilibrium point +X approaches +$. For k = 2 these two points coincide, and the state 
space for n 2 > 2 is depicted in Fig. [|. Both +K_ and +K + still act as past attractors, while _K + is a future attractor. 
However, the stability of has changed; it has become a saddle. There is still a separatrix surface associated with 
+F. Thus, orbits having +K + as their past attractor all end at -K + . 

For k 2 > 37, the equilibrium point +FS, corresponding to the matter-scaling solution, appears from + $, see Fig. 
The point +FS is a spiral sink with an out-going eigenvector direction entering the interior state space. Note that 
this equilibrium point thus is stable in the flat (Q — 1) submanifold, but unstable to curvature perturbations (i.e., 
perturbations in the Qo direction). The scalar-field dominated point still is a saddle, and now there is a separatrix 



i 



G 





7 




FIG. 6. The state space for closed FRW models with a scalar field when k 2 > 6. Dashed curves and white arrows and circles 
are screened. Dotted orbits are in the interior of the state space. 

surface spanned by the out-going eigenvector directions there. Thus, there are two separatrix surfaces, both of which 
are spiraling around the out-going eigenvector direction of + FS. 

When k increases, the + $ equilibrium point comes closer and closer to +K + , and for k 2 = 6 they coincide. The 
state space when k 2 > 6 is given in Fig. |^. There is only one past and one future attractor, namely +K_ and _K + , 
respectively. 

To summarize, when n 2 > 2 all solutions start from and recollape to a singularity (K — > K). Thus, in this case 
solutions can neither expand forever nor inflate. When k 2 < 2, there are also ever-expanding (K — > $) (and ever- 
collapsing $ — > K) solutions in addition to the recollapsing solutions. Inflation occurs when S > 0, i.e., 39 + 9 2 > 0, 
which leads to the condition 



-(37 - 2) - 3(2 - 7 )C/ 2 + 3 7 VF 2 > 0. 



(23) 



This corresponds to a parabolic region along the ridge of the state spaces, Figs. [TJ, || and || - 1]. The only equilibrium 
points within this region are ±$ for k 2 < 2, corresponding to power-law inflation. Consequently, in the case k 2 < 2 
there is a subclass of solutions that inflate. 



III. KANTOWSKI-SACHS MODELS 



We now turn our attention to the Kantowski-Sachs (KS) models. The line element can be written 

ds 2 = -dt 2 + D^tfdx 2 + D 2 (t) 2 dn 2 , 

where 

D 1 = exp [/3°(t) - 2/3+ (t)] , D 2 = exp [f3°(t) + /?+(*)] . 
The kinematic quantities of the fluid congruence are related to the Misner variables (/3°, by 

9 = 3/3°, a+ = 3/3+, 

and the evolution equations for the metric functions B\ = D^ 1 and B 2 = D^ 1 become 



The conservation equations give 



Bx = ~{0 + a+)B u B 2 = --{6 + a+)B 2 . 



/i = -70/i, 



x = -ex + —=v, 

V2 



(24) 
(25) 

(26) 

(27) 

(28) 
(29) 
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and the field equations yield 



2 ^2 



ZB 2 - 3X' 



■SV) 



2o% + 6X Z -W+ -(3 7 — 2)a* J , 



36»cr , 



3X 2 -W- 3/i) 



(30) 
(31) 
(32) 



The Friedmann equation, Eq. (p0|) , together with the assumption /i > shows that D = \/0' 2 + 3B 2 is a dominant 
quantity. Consequently, compact variables are introduced according to 



Qo — jji Q+ — -jj, 



V3X 



W = 



(33) 



D 1 D 

The curvature variable K = 3B 2 9~ 2 = (1 — QqVQo shows that the flat solutions correspond to Qq = 1. The 
Friedmann equation becomes 



n, 



3/i 

TP- 



l-QX-U 2 -W 2 



By introducing a new independent variable, r, where ' = d/dr = d/cft, the evolution equation for D, 



D' 



_(1 - Ql) + 3Q (Q+ + U 2 + ln D } 



D. 



(34) 



(35) 



decouples, and a reduced set of evolution equations is obtained: 



Qd 
Q'+ 



(i-Qo) 



1 + QoQ+ - 3 (Q 2 
-(l-Ql)(l-Q 2 +) + 3Q Q- 

U' = u{(l- Ql)Q+ + 3Q [-1 + (q 2 
W' = W[{1- Ql)Q+ + 3Q (Q+ + £/ 2 



(36) 



-U 2 

In 



-n 



D 



)]} 

3 



kW 2 



kUW. 



There is also an auxiliary evolution equation: 



Q' D = -il D {3 7 Qo - 2 [q+(1 - Qg) + 3Q (O 



[/ 2 



»]} 



(37) 



Note that by setting n = 0, U = 0, and identifying A = Vq, = T'F 2 , the evolution equation corresponding to 
Kantowski-Sachs models with a cosmological constant are obtained Jig. The deceleration parameter is given by 



9pf = ~m 



1-3 Q 



C/ 2 + ^ D 



Note that the dynamical system Eqs. ( |36[ ) is symmetric under the transformation 

(r, Qo, Q+, C/, W) -> (-r, -Qo, -Q+, -£/, ■ 



(38) 



(39) 



Thus, it is sufficient to discuss the behaviour in one part of the state space, the dynamics in the other part being 
obtained by Eq. @. 
The function 



M = Q- 2( ^- 2 \i~Q 2 )- 3 ^n%, 

M' = 2Q; 1 [(3 7 - 2)(1 - Qo) + 3(2 - 7 )Q+] M 
is monotonic in the regions Q + > and Q + < 0, since 2/3 < 7 < 2. Noting that 

Q' + \ Q „ = -(!- Ql)<0, 



(40) 
(41) 

(42) 



we conclude that the submanifold Q + = is not invariant, but acts as a membrane. Thus, the existence of M rules 
out any periodic or recurrent orbits in the interior of the state space and again global results are possible. From the 
expression for the monotonic function M we can immediately see that asymptotically Q+ — » 0, Q 2 , — * 1 or fljj — > 0. 
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Interpretation 


Qo 


0+ 


U 


w 


Note 


±F 


Flat Friedmann 


£ 













±K 


Kinetic dom. 


e 




t/o 





Ring 


= $ 


Scalar-field dom. 
(when re 2 < 2) 




n 






k 2 <r (\ 


±H 


Curvature scaling 


1 i 2 


4+k^ 6 


4+K 2 




< 2 


±FS 


Flat matter scaling 


e 





V 2 « 


\/fiV7(2-7) 


k 2 > 37 


±SSKS 


Self-similar KS 


37-4 6 


3 7 -2 
37-4 c 








7 <2/3 



TABLE V. Equilibrium points of the KS models. 



1 






-37 



-12 



(37-4) 2 





i 

IP 1 + K 2 
1Z (4+ K 2)2 

37 





g P f 

3(37-2) 
|(- 2 -2) 

1 k 2 -2 

2 T+k 2 " 

1(37-2) 
1(37-2) 



TABLE VI. The physical quantities Qd, Q<j>, and q p { for the different equilibrium points of the KS models. 



A. Equilibrium points of the KS dynamical system 

The dynamical system, Eqs. (|36|), has several equilibrium points, which are displayed in Table. As before, e = ±1 
denotes the sign of Qo, while fl^ = U 2 + W 2 . Again, the left subscript gives the sign of Qo and indicates whether the 
corresponding model is expanding or contracting. The values of 17 _p , f^, and q p f for each of the equilibrium points are 
given in Table ^ while the eigenvalues are disp layed in Table VII . Note that all of the equilibrium points correspond 
to exact self-similar cosmological models [p|]lO| . 

As for the closed FRW models, ±F denotes the flat Friedmann solution. Note that the closed FRW solution without 
a scalar field does not appear as a submanifold of the Kantowski-Sachs models without a scalar field. Consequently, 
there is no orbit connecting + F with _F. 

There are two sets ±K of vacuum (fi^ = 0) equilibrium points, parameterized by the constant Uo, corresponding 
to kinetic dominated solutions. These sets are analogues of the "Kasner rings" that are present for various Bianchi 
models. 

The flat scalar-field dominated points ±<I>, already encountered for the closed FRW models, appear in the Kantowski- 
Sachs case as well. As for the closed FRW models, they are physical when n 2 < 6 and inflationary when k 2 < 2. 

There are also equilibrium points ±5, corresponding to curvature scaling solutions (i.e. they have £Id = 0, Q§ < 1) 
which are physical when k 2 < 2. They are also inflationary, but in other respects they resemble the points ±X of the 
closed FRW models. 

As for the closed FRW models, the equilibrium points ±FS, corresponding to the flat matter-scaling solution, enter 
the physical part of the state space when n 2 > 37. 

Finally, there are also equilibrium points corresponding to the self-similar Kantowski-Sachs solution. This solution 
is only physical when 7 < 2/3, and so we will not consider them further. 



The eigenvalues for each of the equilib rium points are given in Table VII. The sources and sinks of the dynamical 



system when 7 > 2/3 are listed in Table VIII (all of the other equilibrium points are saddles). Thus, there is always 
two segments of the equilibrium set + K that act as sources for orbits. Similarly there are two segments on _K that 
are sinks. When k 2 > 2, these are the only attractors, and all solutions start from and recollape to a singularity 
[K -> K). 

When k 2 < 2, which then implies that k 2 < 37 for 7 > 2/3, the equilibrium points ±$ are attractors. Thus, for 
k 2 < 2, there are also ever-expanding ( + K — > + $) and ever-collapsing (_$ — > _K) solutions. 
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Eigenvalues 



±F 


(37 - 2)e 


±K 


4e - 2Q+ 


±<E> 


-I(6- K 2 )e 


±3 


-3f±4e 


±FS 


-§(2-7)£ 


±SSKS 


3 7 -4 6 



-|(2 - 7 )e 
3(2 - 7 )e 
-i(6- K 2 ). 



-6 



7+(7-l)« 
4+k^ 



(37 - 2)e 

37-4 c 



■1(2 -7)e 


-(2- K 2 )e 



3e 



kUq 



-(3 7 -^ 



3 2+ fg 2 ±V(2+^)(l8-7K3) 



& [(2 - 7 )k ± y (2-7)(24 7 2 -(97 -2)^)| e 

3 2-7± % /(2-7)(24 7 2 -41 7 +18) 
2 37-4 6 



TABLE VII. Eigenvalues for the different equilibrium points of the KS models. 



Past attractors 



Expanding from a singularity (Qo > 0) 




+K 


kU < a/6 


Contracting from a dispersed state (Qo < 0) 






k 2 < 2 


Future attractors 


Contracting to a singularity (Qo < 0) 




_K 


kUo < V6 


Expanding to a dispersed state (Qo > 0) 




+$ 


K < 2 



TABLE VIII. Summary of sources and sinks for the KS models. 



From the expression for the monotonic function M we deduce that all orbits asymptotically have Q+ — * 0, Qq — > 1 
or Qjj — > 0. Indeed, the existence of the monotonic function ensures that there are no periodic orbits and that 
generically orbits asymptote towards the local attractors (sinks and sources). Therefore, we can determine the global 
dynamics of the models. 

To summarize, when k 2 > 2 all solutions start from and recollape to a singularity (K — » K). Thus, in this case 
solutions can neither isotropize nor inflate. When k 2 < 2, there are also ever-expanding (K — » $) (and ever-collapsing 
$ — * K) solutions in addition to the recollapsing solutions. Again, the $ points correspond to power-law inflation 
when k 2 < 2. Consequently, in this case there is a subclass of solutions that isotropize and inflate. 

The global asymptotic dynamics is similar to that in the case of positive-curvature FRW models. However, due 
to the presence of shear, the intermediate or transient dynamics can be quite different. In the Kantowski-Sachs case 
the state space is four-dimensional and so we cannot display the phase portraits graphically (as in the FRW case). 
However, as an illustration we shall present the phase portraits in the three-dimensional fluid vacuum and the massless 
scalar field invariant sets in order to compare intermediate behaviours. 



B. Fluid vacuum 



The fluid vacuum (Qrj =0) is an invariant submanifold, as seen from Eq. (|3jj). Using the Friedmann equation to 
eliminate W , we obtain a three-dimensional dynamical system in (Qo, Q+, U): 



Qo = (i - Qo) [i + QoQ+ - ml + u 2 )] , 

Q+ = -(i - Qo)(i - Ql) - 3QoQ+(i - Ql - u 2 ) 



u' = (i - qDq+u 



-k-ZQ U\ {l-QX-U 2 ). 



(43) 
(44) 

(45) 



From table VI. it is immediately seen that the equilibrium points that are contained in this submanifold are ±K, 
and ±S. The state space is depicted in Figs. [?], || and [| Note that Z = 0, where Z is defined by 



1 



Z = Q, 

is an invariant submanifold, and that both ±$ and ±B are contained in this submanifold. 



2^o 



(46) 
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FIG. 7. The state space for Kantowski-Sachs models with no fluid (Q,o =0), with k 2 < 2. Dashed curves and white arrows 
and circles are screened. Dotted orbits lie in the Z = invariant submanifold. 




FIG. 8. The state space for Kantowski-Sachs models with no fluid (Qd = 0), with 2 < k 2 < 6. Dashed curves and white 
arrows and circles are screened. Dotted orbits lie in the Z = invariant submanifold. 
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FIG. 9. The state space for Kantowski-Sachs models with no fluid (Qd = 0), with k 2 > 6. Dashed curves and white arrows 
and circles are screened. 



C. Massless case 



The massless case corresponds to the invariant submanifold W = 0, which leads to a three-dimensional system in 
(Qo,Q+,V): 



Q'o = -(l-Ql) 



%^-Q»Q + + \{2-i){Ql + u*) 



Q'+ = -(l - Qo)(i - Ql) - o( 2 - 7)QoO+fiD. 



t/' = f7 (l-Qg)Q+--(2-7)Q nr. 



(47) 



From table [v], it is immediately seen that the equilibrium points that are contained in this submanifold are ±K and 
±F. The state space is depicted in Fig. HO. 



IV. DISCUSSION 



We have studied closed cosmological models with a perfect fluid satisfying a linear equation of state with 2/3 < 7 < 2 
and a scalar field with an exponential potential. We have utilized a new set of normalised variables which lead to the 
compactification of state space, enabling us to apply the theory of dynamical systems to determine the qualitative 
properties of the models. In all cases we have been able to find monotonic functions which, together with a local 
analysis of the equilibrium points, enable us to determine the global properties of the models. 

We first studied the closed FRW cosmological models. We found that when n 2 > 2, all solutions start from and 
recollape to a singularity (K — > K). In this case solutions generically do not inflate. When k 2 < 2, solutions can either 
recollapse (K — ► K) or expand forever (K — > $) towards power-law inflation solutions (or collapse forever $ — > K); 
consequently, in this case there is a subclass of solutions that inflate. A number of phase portraits were displayed. 

These results generalise previous qualitative work on positive-curvature FRW models with a scalar field (only) || 
and with a scalar field plus a barotropic perfect fluid |g] in which compactified variables were not utilized, and rigorous 
analyses of perfect fluid (only) models using compactified variables []rO| , [l~L| , and completes and generalises more recent 
work using different compactified variables p7| . We also note that positive-curvature FRW models with a perfect 
fluid and a positive cosmological constant have been investigated recently using qualitative methods and utilizing 
compactified variables 

In the case of the Kantowski-Sachs models we again found that when k 2 > 2 all solutions start from and recollape 
to a singularity (K — > K) and can consequently neither isotropize nor inflate. When k 2 < 2, there are also ever- 
expanding (K — * <j>) (and ever-collapsing $ — > K) solutions in addition to the recollapsing solutions, where again the 
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FIG. 10. The state space for Kantowski-Sachs models with fluid and massless scalar field, W = 0. Dashed curves and white 
arrows and circles are screened. The horizontal plane is the [7 = invariant submanifold (compare with pl|). 

4> points correspond to the flat FRW power-law inflationary solution. Consequently, in this case there is a subclass 
of solutions that isotropize and inflate. 

The investigation of Kantowski-Sachs models complements the study of Bianchi models || and completes the 
analysis of spatially homogeneous models. Collins |^] studied perfect fluid Kantowski-Sachs models qualitatively 
using expansion-normalised variables (for which the state space was non-compact) and showed that all models start 
at a Big Bang and recollapse to a final "Big Crunch" singularity. This work was generalised recently by Goliath and 
Ellis |lq| in which Kantowski-Sachs models with a perfect fluid and a cosmological constant were investigated using 
qualitative methods and utilizing the compactified variables of Uggla and Zur-Muhlen Jl3[ ; particular attention was 
focussed upon whether the models isotropize, thereby explaining the presently observed near-isotropy of the universe. 
More importantly, Kantowski-Sachs models with a scalar field and an exponential potential, but without barotropic 
matter, have been studied qualitatively p"^ ], although compactified variables were not utilized. 

To conclude an analysis of positive-curvature spatially homogeneous cosmological models with a perfect fluid and a 
scalar field with an exponential potential, Bianchi type IX models would need to be studied. However, such a study 
is beyond the scope of the current paper. For example, Bianchi type IX models are known to have very complicated 
dynamics, exhibiting the characteristics of chaos |lC|Jl5[l , However, partial results are known. Bianchi type IX models 
with a scalar field (only) have been studied qualitatively, with an emphasis on whether these models can isotropize 
fuj . Scalar-field models with matter have also been studied ||. For example, it has been shown that the power-law 
inflationary solution is an attractor for all initially expanding Bianchi type IX models except for a subclass of the 
models which recollapse ||]4|. However, compact variables have not been utilized and the analyses were not rigorous. 

A more rigorous treatment of the class of Bianchi type IX models with a non-tilted perfect fluid (only) using 
compactified variables has been possible [[lOj . Although an appropriately denned normalised Hubble variable is found 
to be monotonic, enabling some results to be obtained, several problems remain open. More rigorous global results are 
possible. For example, Bianchi type IX models with matter have been shown to obey the "closed universe recollapse" 
conjecture ]23| ], whereby initially expanding models enter a contracting phase and recollapse to a future "Big Crunch". 
In addition, Ringstrom has proven that a curvature invariant is unbounded in the incomplete directions of inextendible 
null geodesies for generic vacuum Bianchi models |24| ] , and rigorously shown that the Mixmaster attractor is the past 
attractor of Bianchi type IX models with an orthogonal perfect fluid (25). A complete qualitative analysis of the 
special class of locally rotationally symmetric Bianchi type IX perfect fluid models, which do not exhibit oscillatory or 
chaotic behaviour near to the initial or final singularities, has been given in based upon an appropriately defined 
set of bounded variables. 

The Kantowski-Sachs models exhibit similar global properties to the positive-curvature FRW models; in particular, 
for k 2 > 2 all initially expanding models reach a maximum expansion and thereafter recollapse, whereas for k 2 < 2 
models generically recollapse or expand forever towards a flat isotropic power-law inflationary solution. The Bianchi 
type IX models share these qualitative properties. However, the intermediate behaviour of the Kantowski-Sachs 
models can be quite different to that of the FRW models. In order to illustrate the possible intermediate dynamics of 
the Kantowski-Sachs models, we studied the special cases of no barotropic fluid, and a massless scalar field in Sees. 
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Ill B and [II C| , respectively (see Figs. fz| — P~0|) - 

Finally, we remark that the dynamics of the Kantowski-Sachs models obtained here will be crucial for understanding 
the global dynamics of general self-similar spherically symmetric models po| . 
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